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The opportunity to study physics at supra-nuclear densities through X-ray observations of neutron 
stars has led to in-depth investigations of certain approximately universal relations that can remove 
degeneracies in pulse profile models. One such set of relations determines all of the multipole 
moments of a neutron star just from the first three (the mass monopole, the current dipole and the 
mass quadrupole moment) approximately independently of the equation of state. These three-hair 
relations were found to hold in neutron stars that rotate rigidly, as is the case in old pulsars, but 
neutron stars can also rotate differentially, as is the case for proto-neutron stars and hypermassive 
transient remnants of binary mergers. We here extend the three-hair relations to differentially 
rotating stars for the first time with a generic rotation law using two approximations: a weak-field 
scheme (an expansion in powers of the neutron star compactness) and a perturbative differential 
rotation scheme (an expansion about rigid rotation). These approximations allow us to analytically 
derive approximately universal relations that allow us to determine all of the multipole moments of 
a (perturbative) differentially rotating star in terms of only the first four moments. These new four- 
hair relations for differentially rotating neutron stars are found to be approximately independent of 
the equation of state to a higher degree than the three-hair relations for uniformly rotating stars. 
Our results can be instrumental in the development of four-hair relations for rapidly differentially 
rotating stars in full General Relativity using numerical simulations. 


I. INTRODUCTION 

The extreme densities inside neutron stars (NSs) have 
sparked interest among many since new physics may lay 
hidden at supra-nuclear densities [1,2]. One tool to probe 
the internal structure of NSs is electromagnetic observa¬ 
tions of binary and isolated pulsars, although gravita¬ 
tional wave (GW) [3-12] and neutrino [2] observations 
may soon be added to the toolbox. X-ray pulse and 
atomic line profiles emitted by millisecond pulsars must 
be fit to models that depend, in part, on intrinsic NS pa¬ 
rameters, such as their mass, radius, moment of inertia 
and quadrupole moment [13-16]. Since the NS equation 
of state (EoS) - the relation between state variables such 
as pressure and density - is unknown, one must fit all 
of these parameters to the data independently , thus di¬ 
luting the information that can be extracted. Reducing 
the number of independent parameters in an EoS inde¬ 
pendent way would help extract more information from 
observables [17, 18], and in particular, it may allow for 
the precise extraction of the NS mass and radius with 
NICER [19] and LOFT [20-22], 

One can partially break the degeneracies among pa¬ 
rameters in NS observations by using universal rela¬ 
tions between certain observables that do not depend 
strongly on the EoS. For example, the moment of in¬ 
ertia (I), the tidal Love number and the quadrupole 
moment (Q or M 2 ) have recently been found to be in¬ 
terrelated in an approximately EoS-independent fash¬ 
ion [23, 24]. The I-Love-Q relations can thus be used to 
reduce the number of model parameters without knowl¬ 
edge of the high-density EoS. Such relations have al¬ 
ready been extended to binary systems with dynamical 
tides [25], proto-neutron stars (PNSs) [26], magnetized 


NSs [2 ], NSs with different EoSs [2 ], rapidly rotating 
NSs [29-33], NSs in a post-Minkowskian expansion [34], 
NSs with anisotropic pressure [35], as well as relations 
with different normalizations [36] , for exotic compact ob¬ 
jects [37] and in modified gravity theories [23, 24, 38-42]. 

The I-Love-Q discovery has also led to the NS three- 
hair relations, which resemble the well-known, black hole 
no-hair relations [43-48] . By using a multipole expansion 
to describe the mass distribution inside a star, its exterior 
gravitational field becomes a function of the mass and 
mass-current multipole moments. Although the expan¬ 
sion involves an infinite sum of moments, the three-hair 
relations reduce all of the moments to functions of only 
the first three in a way that is approximately independent 
of the EoS. Therefore, the NS three-hair relations allow 
one to approximately describe the exterior gravitational 
field of a NS in terms only of its first three, non-vanishing 
multipole moments. 

Such approximate universality is highly sensitive to 
the elliptical isodensity approximation used to derive the 
three-hair relations [49]. This approximation assumes 
that constant density contours of NSs are self-similar el¬ 
lipsoids, with the ellipsoidal radius identical to the ra¬ 
dius of a spherically symmetric star of the same volume. 
The three-hair relations were first derived for weakly- 
gravitating stars [33, 50], and were soon after confirmed 
in full General Relativity (GR) up to hexadecapole or¬ 
der [32], but always assuming uniform (rigid) rotation, 
i.e. assuming a constant angular velocity of rotation. 
This assumption is well-justified for old and cold NSs, 
such as millisecond pulsars in which differential rotation 
has damped out, and it ensures the elliptical isodensity 
approximation holds to a sufficiently good accuracy. 

But there are other physical scenarios in which uniform 
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FIG. 1. (color online). (Left) 3D plot of the four-hair relation showing the normalized mass hexadecapole moment M 4 as a 
function of the normalized mass quadrupole (M 2 ) and mass-current octopole (S 3 ) moments computed with a polytropic EoS 
with various polytropic indices n in the slow-rotation limit. The light blue-green plane represents the fiducial n = 0.65 case. 
Observe that the n = 0.3 (green dots), 0.5 (blue dots), 0.8 (orange dots), 1 (red dots) results are all almost on the same plane. 
The white dashed line in the left panel represents the uniform rotation curve for an n = 1 polytrope. (Right) Contour and 
color gradient plot in the S 3 -M 2 plane of the maximum fractional difference between M 4 computed with a polytropic EoS with 
index n £ [0.3,1] and a fiducial index of n = 0.65 assuming |-y| < 1/2. The red dashed lines demarcate the region of I 7 I < 0.1 
with n= 1. Observe that the maximum fractional difference is always less than ~ 6 % and 3% for all values of S 3 and M 2 
within the |'y| < 1/2 and | 7 | < 0.1 regions respectively. 


rotation is not a good assumption. One such scenario 
are PNSs born after a supernova explosion, which may 
eventually be probed via neutrino and electromagnetic 
observations. Universal relations in PNSs were studied 
in [26], where the original I-Love-Q relations were found 
to hold after only a few hundred milliseconds from birth. 
Such an analysis, however, was limited to uniformly ro¬ 
tating PNSs. Another scenario in which uniform rotation 
is not a good assumption is a hypermassive NS (HMNSs) 
produced after the merger of a NS binary, a prime candi¬ 
date for GW detection [51-5 ]. Differential rotation may 
here be key by providing a temporary mechanism that 
prevents the HMNS from further collapsing. A third sce¬ 
nario in which differential rotation may arise is in NS 
r-modes - toroidal oscillations with a Coriolis restoring 
force - that may be driven unstable via GW emission, 
thus leaving an imprint in GW signals [54], Stars with 
differential rotation have also been found to be more sus¬ 
ceptible to non-axisymmetric instabilities, which can also 
lead to GW signals [55, 56]. 

How does differential rotation affect the three-hair re¬ 
lations? This is the topic of this paper. In order to under¬ 
stand physically the effect of differential rotation in the 
three-hair relations, we work semi-analytically by mak¬ 
ing use of two approximations: a weak-field scheme and 
a perturbative differential rotation scheme. The former 
is simply an expansion in the NS compactness (the ratio 


of its mass to its radius), while the latter is an expansion 
about rigid rotation. For sufficiently small deformations, 
these approximations allow us to continue to use the el¬ 
liptical isodensity approximation; this, in turn, allows us 
to develop an analytic understanding of how differential 
rotation impacts the three-hair relations, as well as to de¬ 
velop new, four-hair relations that are as EoS universal 
as the three-hair relations for rigidly rotating stars. 

Several physical scenarios exist in which the small dif¬ 
ferential rotation approximation is a good approximation 
to Nature, such as during a certain phase after the birth 
of PNSs. When a PNS forms after gravitational collapse, 
it is hot, unstable and highly differentially rotating. Af¬ 
ter tens to a few hundreds of milliseconds, the PNS goes 
into a “quasi-stationary” phase [26], where it is quieter, 
evolves more slowly and can be described as a sequence 
of equilibrium configurations [57-60]. Such a phase lasts 
for roughly one minute, until finally, the star becomes 
cold and uniformly rotating. During this evolution, the 
amount of differential rotation is reduced due to dissi¬ 
pative mechanisms, such as magnetic braking and vis¬ 
cous damping [51, 61, 62]. Therefore, as a PNS evolves 
into a cold NS, an intermediate phase exists (the quasi- 
stationary phase), where the small differential rotation 
approximation is valid. This may not be the case for 
differentially rotating HMNSs formed after binary merg¬ 
ers, which can last ~ 20 ms or longer depending on the 
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EoS, because these can collapse to a BH before reaching 
a phase where the small differential rotation approxima¬ 
tion is valid. 

Even within the small differential rotation approxima¬ 
tion we must still make a choice that determines precisely 
how the star rotates differentially: we need to specify 
a rotation law. We here work with a generic analytic 
parametrization [63] that allows us to model j-constant 
and v-constant laws, the Keplerian angular velocity pro¬ 
file and the rotation law for HMNSs. We further in¬ 
troduce a dimensionless parameter 7 that quantifies the 
extent to which the rotation law differs from a constant, 
with 7 = 0 corresponding to uniform rotation. The rota¬ 
tion law affects directly only the mass-current moments 
Si, because to leading-order in the weak-field expansion 
the mass moments Mi do not depend explicitly on rota¬ 
tion. Differential rotation, however, does affect the stellar 
shape, which does also modify Mi. 

With this at hand, we rederive the three-hair relations 
of [31 ] using a polytropic EoS with polytropic index n. 
We first derive the modified equations of structure ana¬ 
lytically and then solve them numerically for n £ [0.3,1] 
and analytically for n = 0 and for n = e <C 1 in the slow- 
rotation limit, extending Ref. [50] to slowly, differentially 
rotating stars. We then compute various multipole mo¬ 
ments and their interrelations and find that the three- 
hair relations acquire a correction that is proportional to 
7 . This correction is mildly EoS dependent, deteriorat¬ 
ing the universality of the three-hair relations by at most 
0 {7 x 10%). But since I 7 I -C 1, this deterioration is al¬ 
ways much less than 10 %, and thus, much smaller than 
the inherent EoS variability of the three-hair relations in 
the rigidly rotating case. 

These semi-analytic results suggest the construction 
of new four-hair relations for differentially rotating NSs 
by replacing the 7 dependence of the three-hair rela¬ 
tions with the next, independent multipole moment (the 
mass-current octopole S3). The resulting four-hair rela¬ 
tions are then independent of 7 and EoS-independent to 
approximately 0 ( 10 %), just as in the three-hair case of 
rigidly rotating stars. The left panel of Fig. 1 presents 
an example of the four-hair relations, where we plot M 4 
in terms of M2 and S3, normalized to their correspond¬ 
ing black hole expressions and using a polytropic EoS 
with n £ [0.3,1]. Observe that regardless the value of 
n , the multipole moments lie approximately on the same 
so-called invariant, four-hair plane. The right panel of 
Fig. 1 shows the maximum fractional difference in the 
S3-M2 plane between M 4 computed with a polytropic 
EoS with n £ [0.3,1] and a fiducial EoS with n = 0.65, 
assuming [yj < 1/2. Observe that the maximum frac¬ 
tional difference stays always below ~ 6 %. 

How does the EoS variation in the four-hair relations 
compare to the three-hair relations for uniformly rotating 
stars? The I 7 I < 0.1 region of the right panel of Fig. 1 
shows that the maximum EoS variation in the four-hair 
relations is ~ 3%, which is actually smaller than the 
~ 4% variation of the three-hair relations for uniformly 


rotating stars [33]. This means that small amounts of 
differential rotation actually improve the universality, al¬ 
beit by a small amount. The reason for this improvement 
is that the four-hair relations contain an additional de¬ 
gree of freedom (S3) relative to the three-hair relations. 
The EoS variation can then be explored along a new di¬ 
rection (with a fixed M 2 and S 3 ) that is not allowed in 
the three hair case and that can be exploited to minimize 
the EoS variation further. Of course, for large differen¬ 
tial rotation, i.e. I7I > 0.1, the EoS variation increases, 
and eventually, the EoS variation in the four-hair relation 
exceeds that in the three-hair relations with uniform ro¬ 
tation. 

The remainder of this paper presents the computa¬ 
tional details that lead to the above results. Section II 
introduces the generalized differential rotation law we use 
in this paper and the concept of multipole moments for 
NSs. Section III derives the three- and four-hair relations 
of differentially rotating stars and discusses their EoS- 
dependence. Section IV concludes and points to future 
research. Henceforth, we follow the notation of [33], and 
in particular, we use geometric units in which G = 1 = c. 

II. MULTIPOLE MOMENTS FOR 
DIFFERENTIALLY ROTATING STARS 

The gravitational potential around a source, assuming 
an asymptotically flat and axisymmetric spacetime, can 
be completely described by its mass and mass-current 
moments. To leading- (so-called Newtonian) order in a 
weak-field expansion (C := M*/R* <C 1, with M* the 
mass and I?* the equatorial radius), the multipole mo¬ 
ments of differentially rotating stars are given by [64] 

r 1 rR*( m) 

Mt = 27 t J J p(r,p)r e+2 dr P e (p)dp, (1) 
Si = J ^ fl{r,p)p{r,p)r e+3 dr 

X ^M(1 -p*) d p, (2) 

dp. 

where p = cos 6 , R*(p) is the stellar surface profile, 
p{r,p) is the star’s density, Pi{p) are Legendre poly¬ 
nomials, and tt(r,p) is the angular velocity of the star. 
Clearly, when Q(r,p) = fl c = const., the star rotates 
uniformly. In the weak-field limit, the gravitational po¬ 
tential of axially symmetric stars can be completely de¬ 
scribed by the mass multipole moments. In general, re¬ 
flection symmetry about the stellar equator forces half of 
the moments to vanish: M 2 ^ + 1 = 0 = S 2 ^. 

Differential rotation modifies the multipole moments 
in two different ways. First, the radial and polar an¬ 
gle dependence of O directly modifies the mass-current 
moments. This dependence is encoded in the so-called 
rotation-law, for which we here adopt the generalized 
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model of [63]: 


^ = ( 1 - a7 ^ si,,a ")’ • (3) 

where fl c is the angular velocity at the core and a\ is 
the semi-major axis of the uniform rotation configura¬ 
tion. The parameter a determines the type of rotation 
law. For the j-constant and v-constant laws, the Ke- 
plerian angular velocity profile and the rotation law for 
HMNSs, a = —1, —2, —4/3 and —4 respectively [63]. The 
dimensionless parameter 7 controls the amount of dif¬ 
ferential rotation, with uniform rotation recovered when 
7 = 0, and with a and 7 satisfying ay < 0 [ ]. The 
perturbative differential rotation approximation consists 
of assuming that |y| <C 1 and expanding all equations in 
this quantity. Such an expansion is a perturbation about 
the uniform rotation case, reducing the rotation law to 

^ = l-7^sin 2 0 + 0( 7 2 ). (4) 

ll c ai z 

Notice that the a has canceled out. 

The second modification that differential rotation in¬ 
troduces is a change in the stellar shape i?*(/r). Fol¬ 
lowing [33], we assume that 7?*(/r) corresponds to the 
radial profile of a constant density star. In the uniformly 
rotating case, this corresponds to spheroids, but when 
when differential rotation is present the shape is no longer 
spheroidal. The change in stellar shape can be expressed 
as an infinite sum of even Legendre polynomials Pe(-)'. 

R*(fJ-) = a 1 ^1+ j^j^V+7 X^ P 2 j(M)a 2 i(e 2 )j 

+ 0{ 7 2 ), (5) 

where a 2 j{e 2 ) is a function of the stellar eccentricity 
squared, e 2 = 1 — a§/af, with 03 the semi-minor axis 
of the uniformly rotating configuration. Appendix A 
presents a mathematical derivation of this result. 

Further assumptions and approximations allow us to 
make more analytical progress. Following [65], we adopt 
the elliptical isodensity approximation, where we assume 
the layers of constant density are ellipsoidal with the 
same eccentricity throughout the star. In reality, iso¬ 
density layers become more spherical towards the center, 
as illustrated in Figure 2, but this has been shown to 
have a very small effect in the structure of uniformly 
rotating NSs, i.e. calculations that employ the elliptical 
isodensity approximation reproduce full numerical calcu¬ 
lations up to the mass-shedding limit with errors of at 
most 3% [65]. The elliptical isodensity approximation 
introduces an even smaller error when computing multi¬ 
pole moments, since the mass near the center of a star 
contributes less than the mass at larger radii [24]. 

Although the elliptical isodensity approximation has 
only been shown to be accurate for uniformly rotating 
stars, one can still impose such an approximation for 


differentially rotating stars as long as the amount of 
differential rotation is small. In order to explain this 
point further, let us introduce a dimensionless parame¬ 
ter e (<C 1 ) that characterizes the difference between the 
three-hair relations among multipole moments for uni¬ 
formly rotating stars with and without imposing the ap¬ 
proximation. Then, expanding the three-hair relations 
for differentially rotating stars around both 7 = 0 and 
e = 0 , one can decompose the three-hair relations or¬ 
der by order, obtaining a bivariate series with structure 
e>( 7 °e°) + O^e 1 ) +O(~i 1 e 0 ) + 0{j 1 e 1 ) + 0{^ 2 , e 2 ). The 
first term is the three-hair relations for uniformly rotat¬ 
ing stars within the elliptical isodensity approximation. 
The second term is the correction to this due to the 
breakage of the approximation, which has been shown 
to be small [32]. The third term is the three-hair rela¬ 
tions for differentially rotating stars within the ellipti¬ 
cal isodensity approximation, which we calculate in this 
paper. The fourth term is the correction to the latter 
due to the breakage of the approximation. But notice 
that this fourth term is of higher order than O( 7 1 e 0 ) 
by O(e). Therefore, we can safely impose the elliptical 
isodensity approximation even for differentially rotating 
stars as long as differential rotation is small and we per¬ 
form all calculations perturbatively. 
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FIG. 2. (color online). Stellar density (color gradient) and 
isodensity layers under uniform rotation in the elliptical iso¬ 
density approximation (solid contours) and for a typical NS 
(dashed contours) constructed with a realistic SLy EoS [66]. 
The color gradient legend corresponds to the ratio of ellipti¬ 
cal radius f to the stellar surface 7?*. Observe how the solid 
contours are close to the dashed contours, except close to the 
stellar core. 


The elliptical isodensity approximation suggests that 
to solve Eqs. (1) and (2) we should introduce an 
adapted coordinate system x 1 = f 0 (cos 0 ) n l where n l = 
(sin 9 cos </, sin 9 sin <j >, cos 9) is the unit direction vector in 
spherical coordinates. The change in stellar shape man- 
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ifests as a correction to 0(/z), namely 
I 2 00 

©( m ) = \J 1 + +'Yj2 F>2 i(t J ')a2j{ e2 ) + 0(l 2 ), 

j —0 

( 6 ) 

where we used Eq. (5) and the stellar surface is defined 
at r = a\. 

Changing coordinates allows the multipole moment in¬ 
tegrals to be separated into radial and angular parts. The 
mass multipole moments, Me, become 


work could easily extend these results to more compli¬ 
cated EoSs [2, 68]. We first transform to dimensionless 
variables, defined via p = p c [$(£)]"> where £ = (£i/ai)r 
is a dimensionless radius, such that £ = £1 corresponds 
to the stellar surface, and p c is the central density. 
After this transformation, the radial integral becomes 
Ri = Pc(ai/€iY +3 R n ,e, where [33] 

n n ,i= (15) 

J 0 


M e 


2irRe 


h,3 +7(7 + 3)7^2 


+ 0 { 7 2 ), 


(7) 


where [33] 

rai r 1 

Re= p(r)f e+2 df , h,k = / Pe{p)Q(p) e+k dp. 
Jo J -1 

( 8 ) 

Notice the density depends only on f, and the deviation 
from spheroidicity is encoded in It , 2 as 

/ 1 00 

P e (p)Q(p,y +k P2j(p)a 2j (e 2 )dp . (9) 

- 1 3 =0 


Using Eq. (4) to model differential rotation, the mass- 
current moments, Si, become 


Se 


27+1 ' 
7(7 +1) 


a?(27+l) 


Re+i {si e , 3 +7 

(Sit, 5 — Sle^s) 


(7 + 4)57^2 

lO( 7 2 ), (10) 


where 6 N l>k := Nt- 1++2 - Nt+i,k for N = 7, 7, and 7 
with the last defined via 


7^,fc = / p 2 Pe(p)Q(pY +k dp 


( 11 ) 


/-1 


The closed-form solutions to It,k and It,k for various k 
values can be derived using Eq. (7.226.1) from [67]; from 
these, we find 


Sit, 3 = (-1) 

Sit, 5 = (-1) 
x[(e 2 

Sit,. 5 = (-1) 

x[(e 2 


^ 2(2 £+ 1) wr^ e ^ 

7(7 + 2) V ’ 

7(7 + 2)(7 + 4) V ' 
-2)7-3] , 

f±3 2(27 + 1) / _ 2 \ 3 / 2 pt —3 

7(7 + 2)(7 + 4) ' 

-1)7 + 1] , 


( 12 ) 


(13) 


(14) 


while SIe,k will be solved later under the slow-rotation 
approximation. 

By separating the integrals, all of the EoS-dependence 
that remains is in the radial integral, Re. Let us then 
compute these quantities using a polytropic EoS, p = 
Kp 1+l l n , where p is pressure and p is density; future 


The spherically-symmetric, Lane-Emden function $ sp h 
has replaced d under the assumptions of the elliptical 
isodensity approximation, i.e. the constant f surfaces are 
assumed to have the same density profile as a spherically 
symmetric star of the same mass. We drop the subscript 
“sph” below for simplicity. 

Now that the radial integrals have been normalized for 
a polytropic EoS, the multipole moments become 


Me = M 0 


<A h,3 Rn,t 

epTo 7 3 WM 

Oil 2 ) , 


1 + 7 ( (7 + 3) — 3-Y— 

\ 1 t,3 1 0,3 . 


(16) 


27MoO c aq +1 Sit, 2 , !Z n ,t+i 

(27+1)^+ 3 7o, 3 |0'(£i)I 


S t = 


1+7 


(7 + 4) 


Sit ,2 
Sit ,3 


- 3 


To,2 7(7 + 1) (Sit ,5 — SR, 5 ) IZ n ,e + 3 
I 0 ,3 (27+1) £? SR ,3 Hn,t+ 1 


0( 7 2 ), 
(17) 


where one can solve for It ,3 to find 

j4 -3 = (-i) l /TT v ^ e *- (18) 


In deriving these expressions we have used Eq. (7) with 
7 = 0 to solve for p c in terms of the mass monopole Mo, 
where |'d , (£i)| is the derivative of a dimensionless func¬ 
tion related to density at the stellar surface. Notice the 
differential rotation term for the mass moments does not 
depend on the EoS, unlike the mass-current moments. 


III. UNIVERSAL RELATIONS FOR 
DIFFERENTIALLY ROTATING STARS 

A. Three-Hair and Four-Hair Relations 


Let us begin deriving the universal three- and four-hair 
relations by normalizing the multipole moments, using 
the same normalization of [33]: 


Mt = (-) = 


Me 

M^+V ’ 


St = (—)"®~ 


5, 

M<+V ’ 


(19) 


where \ — S\/Mq. It should be noted that M 0 = S\ = 1 
since Mq and S± are used to normalize the moments. The 
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normalized moments M 2 e +2 and S 2 f_+i are given explic¬ 
itly in Eqs. (Bl) and (B2) of Appendix B respectively. 
Inspecting these equations, one notices that all of the 
EoS-dependence in the differential rotation corrections 
can be absorbed into a coefficient, 


C n ,e '■= 


Hn, 2^+4 
£1 Tin, 21+2 


( 20 ) 


Therefore, by studying how much C n / varies with n for 
a set of t, we can determine how much that mass and 
mass-current multipole moments vary with the EoS. 

Deriving the first universal relation begins by solving 
M 2 for a\ and substituting this into Mit+i/Sit+i- Re¬ 
taining only terms up to 0 ( 7 ), one finds 

M 2 e +2 = M 2 S 2 (.. 1-1 (l + 7 Q-n/') + ^(7 2 ) ! (21) 


where the differential rotation term ap £ depends on C n e 
and can be found in Eq. (B3). Notice that Eq. (21) agrees 
with the uniform rotation result of [33] when 7 = 0 . 

Deriving the second universal relation begins by solv¬ 
ing S 2 e+i for a\ l and inserting af f (»§ 2 £+i) into M 2 t+ 2 - 
Retaining terms up to 0 ( 7 ), one finds 


M 2e+2 = A n j(S 2 t+i) 1+1/e ( 1+7 «S) + 0(7 2 ), (22) 

where a second coefficient defined in [ 

_ _(2£+s (K^k-%+2) 

An ’ ( = 3 1 + 1 A£2|tf/(£ 1 )| (23) 

absorbs the EoS-dependence. The differential rotation 
term cv t can be found in Eq. (B4) . 

By combining these universal relations, all of the multi¬ 
pole moments are related through the three-hair relations 


M< 


2C+2 ' 


+ iS 2 ( + i — B n jM. 


n,l ivl 2 


M 2 (l + 7 af’ M) j 


+iS 1 (l+'yatff)] +0( 7 2 ), (24) 


where a third coefficient, also from [33], 



‘ U.O \J.\J 

n n 


FIG. 3. (color online). B n ^ and C n ,e are plotted above for 
the first four £’s with n £ [0.3,1]. The solid lines represent 
the exact numerical results, while the dashed lines come from 
the perturbative approach. The bottom panels show the frac¬ 
tional difference of the semi-analytic results using a fiducial 
n = 0.65. 


where we have defined the piecewise function 


a 


( 3 ) 



even£ 


odd£ 


and where iq = \JM 2 /M 0 , a = S\/M 0 and |_-J is the 
floor function. The three-hair relations closely resemble 
the well-known no-hair relations for black holes: Mf H + 
iSf a = M 0 (i a y [47]. 

Replacing the differential rotation parameter 7 with 
the next independent multipole moment S 3 turns the 
three-hair relations into a four-hair relation. By taking 
the imaginary component of Eq. (24) with i = 1, the ex¬ 
pression can be solved for 7 in terms of M 2 and S 3 . The 
result simplifies into relations that only depend on the 
first four multipole moments: 


B n ,v 


3*+ 1 £ 1 2 V (£ 1 )|^ n , 2£+2 

2f + 3 


(25) 


again absorbs the EoS-dependence. The differential rota¬ 
tion corrections of ’ M ' > and ot~pp are given in Eqs. (B5) 

and (B 6 ) respectively. Notice that ap M ' > does not de- 
(3 S) 

pend on the EoS while or n does. Notice also that these 
relations match those from [33], but also include correc¬ 
tions due to differential rotation to leading-order in 7 . 

For comparison, we temporarily restore units in the 
three-hair relations: 


M 2 f +2 + iS 2 t+i — Bn^M^ M 2 ^1 + cp’P^ 


+iSi (l + «i/ S) ) 


n,£ 

oh 2 ), 


(27) 


where both PpP ' 1 and ot~pp depend on M 2 and S 3 (but 
not on 7 ), and the moments were normalized by Mq and 
S±. These differential rotation corrections can be found 
in Eqs. (B7) and (B 8 ) of Appendix B. 


B. Generic EoS Dependence 


Me + i-Se = B n ^e~i^AI n 


(iqf (l + l a n/) + 0 ( 7 2 ), 


( 26 ) 


We have found differential rotation-corrected three- 
hair relations and new four hair relations, but how EoS 
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independent are these? Part of the EoS-dependence is 
absorbed in the three coefficients A n /, B n /, and C n /, so 
we can analyze the variation of these with the polytropic 
index n. We will only look at B n / and C n /, since the 
EoS-dependence of A n j was already analyzed in [33] in 
the uniformly rotating case, and found to lead to EoS 
variability of less than 10%. 

Figure 3 shows the relative fractional differences of 
B n / and C n / computed with a polytropic EoS with 
n € [0.3,1] and a fiducial polytrope with n = 0.65. As 
found in [33], observe that the EoS variation in the uni¬ 
form rotation coefficient, B n /, increases as one increases 
£. However, B ni 3 corresponds to the relation between 
Sj or M$ against M2 • Therefore the ~ 40% variation in 
B„ 3 between n = 0.3 and 1 are for less significant, higher 
order moments. Observe also that the EoS-dependence 
of the differential rotation correction C n / is smaller than 
that of the uniform rotation term B n / for large £, while 
the opposite is true for small £. However, recall from 
Eq. (27) that C n / is multiplied by 7, which is assumed 
to be much less than unity. Therefore, differential rota¬ 
tion does not destroy the EoS universality of the three- 
hair relation among lower multipole moments, at least 
when treated perturbatively and to Newtonian order in 
a weak-field expansion. 

We can gain some further mathematical understand¬ 
ing of this result by considering perturbations about an 
n = 0 polytrope, for which we can solve the Lane-Emden 
(LE) equation exactly. Let us then set n = h + e, where 
h is a background polytropic index for which we assume 
we have a solution to the LE equation, and e is a pertur¬ 
bation to n. The solution to the LE equation must then 
become •& = 7) + e57) + 0(e 2 ), where 7) is the background 
solution to the LE equation and 57) is a perturbation. By 
substituting the perturbed variables into the LE equation 
and expanding to leading order in e, the perturbed LE 
equation becomes [50, 69-71] 


With all of these analytical results at hand, we can 
now evaluate the dependence of the EoS factors B n / and 
C n / as functions of n. First, we expand Bo/ = Bo/ + 

eSBo/ + 0 (e 2 ) and Co/ = Co/ + eSCo/ + 0 (e 2 ), noting 
that the background coefficients are given by 


B 0 ,e = ^ 


15 

(2£ + 3)(2£ + 5) J ’ 


A 21 + 5 

o,e 21 + 7' 


(31) 


The perturbation can then be solved for to find 


5B 0 / 

Bo/ 

SCp/ 

Co/ 


-1 

15 


15 H ( £ + - ) - 6 £ - 46 + 15 ln(4) 


(32) 


-2 

21+7 ’ 


(33) 


where H{£) = 1 /k is the £th harmonic number. 

The same derivation for 5Aq//A o/ was presented in [50], 
where the only change from Eq. (32) is a factor of — l/£. 
The EoS-dependence of the coefficients can be analyzed 
by varying e. 

Figure 3 compares these perturbative analytic solu¬ 
tions to the numerically calculated B n / and C n / in the 
top panels. The largest fractional difference between the 
numerical and perturbative calculations occur for B\ o 
(a 10% error) and for C\ t o (a 2.6% error). We then see 
that the perturbative calculation is more accurate for the 
differential rotation coefficient than for the uniform rota¬ 
tion one. By comparing Eqs. (32) and (33) we can also 
see how 5B 0 / is typically larger than SC 0 / for large £, as 
observed in Figure 3. 


C. EoS Dependence in the Slow-Rotation Limit 


^ + ln(i?)^ +1 e = 0 . 


(28) 


We require ffi?(0) = 0 = <5$'(0), because 1 ) already satis¬ 
fies the initial conditions of the LE differential equation. 

The perturbed LE equation can be solved exactly when 
the background polytropic index is n = 0. Doing so, we 
find [50, 69] 
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5 — 3 In 


1 -? 

6 


+ 3 In 1 


4\/6tanh 1 (£/V6) 

£ 


-4. 



From this solution, we can determine the perturbation to 
the location of the NS surface. Letting £1 = £1 + e<5£i + 
0 (e 2 ), where £1 = \/6 and requiring that $(£1) = 0, we 
find the perturbation to the NS surface to be [50, 6! 


(30) 


In the previous subsection, we investigated the EoS- 
dependence of the differential rotation correction C n / of 
Eq. (20), which clearly does not depend on the angular 
integrals le.k of Eq. (9). These angular integrals, how¬ 
ever, enter the three- and four-hair relations in Eqs. (24) 

and (27) through the coefficients. We see that 

they depend on the Legendre coefficients a.ij that charac¬ 
terize the deformations of the stellar surface. Appendix A 
derives such a deformation and shows how to solve for the 
Legendre coefficients in the slow-rotation limit, i.e. ex¬ 
panding in x '■= Si/Mq <C 1. The first three a2j’s are 
given explicitly in Eq. (A24). 

In the slow-rotation approximation, we can then com¬ 
pute the first few multipole moments from Eqs. (16) 
and (17). The first mass moments are 


M 2 

m 4 


ale 2 M 0 K n ,2 ( 16'\ 

" 3£f|i9 , (£i)| V 7 5 ) 

aje 4 MoK nA / _ 
5£^'(6)l V e 2 15/ 


+ 0 ( 7 2 ), 
+ 0 ( 7 2 ), 


(34) 

(35) 


<5^=° = v / 61n(4) 


7 

71' 

















100 



m 2 


FIG. 4. (color online). Contour plot of 7 ( 11 / 2 , 53 ) under 
the slow-rotation limit. Observe how M 2 and S 3 specify the 
amount of differential rotation, and thus, if | 7 | <0.1 then 
only a subset of M 2 and S3 are allowed. The contours are 
generated with a polytropic EoS with n = 1 because values 
of n < 1 yield less restrictive contours and thus a wider range 
of allowed M 2 and S3. The white dashed zero-contour cor¬ 
responds to the uniform rotation case, which shows a linear 
dependence between M 2 and S3 that agrees with the three- 
hair relations for an n = 1 polytrope. The green lines at 
M 2 = 30 and S3 = 80 show the range of S3 and M 2 used in 
Fig. 5 that is kept within | 7 | < 1/2, demarcated by the blue 
lines. 


and the first current moments are 


Si 

S 3 


S 5 


)+0(7 2 ), (36) 

2afe 2 M 0 fl c TZ nA / 7 24 TZ n>e \ n( 2 , 

5£?TO)I 1 e2 49^ft n J + 17 

(37) 

2afe 4 M 0 fl c 7?.„, 6 I" 78/35 77.„, 8 \ 

7£?|0'(£i)l L + 3gKZe). 

0( 7 2 ) • (38) 


2a 2 AI 0 V, c TZ n2 ( ^ 8 7?- rii 4 

3^1 |^'(^l)| 'v " 7 15^2 Tin ,2 


Notice that the differential rotation corrections to M 2 
and Si are proportional to 7 while those to M4, S3 and 
S5 are proportional to 7 /e 2 . This does not mean that 
our expressions are divergent in the e 2 —> 0 limit; if 
one expands out the expressions above, one finds that 
they are all finite in this limit. The 7 /e 2 dependence, 
however, does mean that our double expansion of slow- 
rotation and small differential rotation is only valid when 
I 7 I C e 2 < 1. Notice that S 3 oc il 3 + jfl, and thus, it is 
non-vanishing to linear order in spin in differentially ro¬ 


tating stars, as already implied in Refs. [ 6 , 72-74] for a 
j-constant law in full GR. Notice also that M 4 oc fl 4 +7 fl 2 
and S 5 oc Q 5 + jfl 3 , and thus these moments are non¬ 
zero at second and third order in spin respectively for 
differentially rotating stars. 

From Eqs. (34)-(38), we can further compute the four- 
hair relation in the slow-rotation limit to find 


M 4 = M 2 B nA 


1 49 (S 3 — M 2 B n ^ 1 ) 

+ 45 M 2 B nil C ntl 


+ 0 ( 7 2 ), (39) 


5s = M 2 B n 


o(i 2 ). 


1 + 


49 (9 + 25C TIj2 ) (S 3 — M 2 B n i) 


45 
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M 2 B n ,C n i 


(40) 


Observe how these relations depend on C n j and B n j. as 
discussed in the previous subsection. 

Before we can plot these four-hair relations in the slow- 
rotation limit, we must choose a range of values for M 2 
and S3. In doing so, however, one must be careful to pick 
a range that satisfies jyl < 1 so that the perturbative 
differential rotation approximation is not violated. The 
dependence of M 2 and S3 on 7 is shown through con¬ 
tours in Fig. 4. Observe that only a subset of values of 
( M 2 ,S 3 ), shown with dotted red lines, lead to I 7 I < 0.1. 
This range of values is in fact the smallest range possi¬ 
ble, since the contours were constructed with a polytropic 
EoS with n = 1 and other values of n < 1 would allow 
for a larger range. 

With the range of allowed values explored, we now 
plot the four-hair relations in the slow-rotation limit in 
two dimensions (Fig. 5) and in three-dimensions (Figs. 1 
and 6 ). The left panel of Fig. 5 shows the M 4 -M 2 relation 
and the S 3 -M 2 relation for a fixed value of S3 = 80 over a 
range of values of M 2 that guarantee that I 7 I < 1/2 (the 
horizontal green line in Fig. 4). Similarly, the right panel 
of Fig. 5 shows the M4-S3 relation and the S3-S3 relation 
for a fixed value of M 2 = 30 over a range of values of S3 
that again satisfy jy] < 1/2 (the vertical green line in 
Fig. 4). Observe how the relative fractional difference 
(bottom panels) shows a weak EoS-dependence, i.e. a 
weak dependence on n, relative to a fiducial polytropic 
index of n = 0.65. The maximum fractional differences in 
the M a -M 2 and M 4 -S 3 relations is ~ 7%, while in the S 5 - 
M 2 and S3-S3 relations it is ~ 20%. But notice that these 
differences decrease to ~ 3% and ~ 10% respectively for 
|y| < 0 . 1 , as shown by the horizontal black lines; again, 
differential rotation does not spoil the approximate EoS 
universality. 

Let us now compare the EoS variation in the four-hair 
relations for differentially rotating stars to that in the 
three-hair relations for uniformly rotating ones. The lat¬ 
ter is simple to obtain: we simply take the difference 
between M 4 using n = 0.65 and n = 1 polytropes with a 
fixed value of M 2 , shown by the vertical magenta line be¬ 
tween the magenta squares in the inset of the top, right 
panel of Fig. 5, a zoom to the region around S3 = 71 
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FIG. 5. (color online). (Top panels) Four-hair relations for differentially rotating stars. The normalized moments M 4 and £5 
are plotted against M 2 with £3 = 80 (left) and S 3 with M 2 = 30 (right) for various polytropic indices n in the slow-rotation 
limit. Magenta squares represent the uniform rotation case (7 = 0 ), while the black diamonds show where 7 = 0.1 for a given 
value of n. The fiducial polytropic index of n = 0.65 is plotted as a solid black line. The inset zooms into a region around 
S 3 = 71 and M 4 = 2.1 x 10 3 , where the vertical magenta (black) line corresponds to the difference in Mj between an n = 0.65 
and n = 1 polytrope for uniformly (differentially) rotating stars with fixed JVf 2 - Observe that the latter is smaller than the 
former when | 7 | < 0.1 with an n = 1 polytrope. (Bottom panels) The fractional difference of each of the four-hair relations 
relative to that of a fiducial polytrope (n = 0.65). Observe that the maximum EoS variation is roughly 7% for M 4 and 20% 
for £ 5 , but this goes down to 3% for M 4 and 10% for £5 in the |y| <0.1 region marked by the horizontal black lines. The 
horizontal magenta lines show the maximum fractional difference in the three-hair relations for uniformly rotating stars. 


and M 4 = 2.1 x 10 3 . The EoS variation in the differen¬ 
tially rotating case, however, is more difficult to obtain 
because it depends both on M 2 and S 3 . For example, if 
we fix M 2 , as we do in the top, right panel of Fig. 5, the 
EoS variation between n = 0.65 and n = 1 polytropes is 
represented by the distance between the black solid line 
and the red dashed line, which clearly depends on the 
value of S 3 chosen. Let us fix the value of S 3 such that 
the EoS variation is largest within the range of values of 
S 3 that corresponds to |y| < 0.1 for an n = 1 polytrope, 
i.e. the largest vertical distance between the black solid 
line and the red dashed line within the region delimited 
by the black dots of the red dashed line. 

With this at hand, we can now compare the EoS vari¬ 
ations, shown in the bottom, right panels of Fig. 5. The 
horizontal magenta line shows the relative fractional dif¬ 
ference in the uniform rotation case, while the black solid 
line shows the maximum, relative fractional difference 
in the differential rotation case for |y| < 0.1. Observe 
that the universality improves slightly when one includes 
differential rotation, i.e. the horizontal magenta line is 
slightly above the horizontal black line. This result is 
robust to the choices made above, since if we had picked 
different values of S 3 (within |'y| < 0.1 for an n = 1 poly¬ 
trope) to compute the EoS variation in the differential 
rotation case, the latter would have been even smaller, 
i.e. the horizontal black line would have been farther 
down. The reason for this is that differential rotation 
introduces a new degree of freedom, encoded in S3 in the 


four-hair relations, which allows the EoS variation to ex¬ 
plore a new direction (the vertical direction in the inset 
of Fig. 5). For small values of I 7 I, there are values of 
S 3 , different from those required in uniform rotation for 
a fixed M 2 , that lead to a smaller degree of variation. 
Of course, if one considers |y| > 0.1, the EoS variation 
in the four-hair relations eventually exceeds that in the 
three-hair relations for uniformly rotating stars. Similar 
features can be seen for the S 3 -S 3 case with fixed M 2 . For 
the S 5 -M 2 case with fixed S 3 , differential rotation dete¬ 
riorates the universality from the uniform rotation case 
slightly, but the latter does not correspond to the uni¬ 
versality of the original three-liair relations in [33], since 
that fixes M 2 instead of S 3 . 

Three-dimensional plots of the four-hair relations in 
the slow-rotation limit are shown in Figs. 1 and 6 . The 
left panel of Fig. 6 shows how the S 5 -M 2 -S 3 relation ex¬ 
ists on an invariant plane that is approximately insensi¬ 
tive to the EoS. The right panel shows this insensitivity 
in more detail, by plotting the maximum fractional dif¬ 
ference in S 5 in the part of the S 3 -M 2 plane that satisfies 
[ 7 ! < 1/2. The yellow (green) solid line shows where 
the fractional difference of £5 computed with an n = 0.3 
(n = 1) polytrope and a fiducial n = 0.65 polytrope van¬ 
ishes. Observe how these lines lie inside the purple region, 
where the EoS variation is smallest. Observe again how 
the maximum EoS variability is always less than ~ 20% 
inside the I 7 I < 1/2 region, but it is less than ~ 7% in 
the I 7 I < 0.1 region. 
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FIG. 6. (color online). Same as Fig. 1 but for S$ as a function of M 2 and S 3 . The yellow (green) line in the right panel 
shows where the fractional difference of S 5 computed with an n = 0.3 (n = 1) polytrope vanishes with respect to the fiducial 
n = 0.65 one. Observe that the maximum fractional difference is always less than ~ 20% for all values of S 3 and M 2 that 
satisfy |y| < 1/2, while it is only between 3% and 7% in the |^y| < 0.1 region (dotted red lines). 


IV. DISCUSSION 

We have studied the effect of differential rotation in 
the approximately universal three-hair relations for neu¬ 
tron stars. To allow for a semi-analytic treatment, we 
used a leading-order weak-field, perturbative differential 
rotation and slow-rotation expansion. Under these ap¬ 
proximations, our results show that differential rotation 
affects the three-hair relations only mildly, without spoil¬ 
ing the universality found in uniform rotation for suffi¬ 
ciently weak differential rotation. We also found that 
differential rotation introduces a new parameter into the 
three-hair relations that quantifies the degree of differ¬ 
ential rotation. We eliminated this parameter in favor 
of the next non-vanishing multipole moment (the mass- 
current octopole) to derive new four-hair relations, which 
typically have smaller EoS variation for small differential 
rotation than the three-hair relations for uniformly ro¬ 
tating stars. 

In order to proceed with our analysis, we chose a 
generic class of rotation laws that encompasses the j- 
and v- constant laws, as well as Keplerian and HMNS 
rotation laws, all of which are identical to linear order 
in small differential rotation. Although other classes of 
rotation laws could have been used to model more re¬ 
alistic differential rotation curves, we expect such laws 
to only induce corrections to the three-hair relations of 
0 ( 7 ) relative to the uniformly rotating case, as long as 
one imposes the small differential rotation approxima¬ 
tion. Therefore, provided that the amount of differential 
rotation is small, the results presented here should also 


be valid qualitatively for other classes of rotation laws. 

Our results suggest that the approximate EoS univer¬ 
sality may survive the presence of differential rotation, 
but a much more detailed numerical analysis would be re¬ 
quired to confirm this. In particular, one could consider 
the multipole moments in full GR (without a weak-field 
expansion), for rapidly rotating stars (without assuming 
slow-rotation) and for large differential rotation, such as 
that produced in the remnant shortly after a supernova 
or a binary merger. Doing so, one could extend the work 
of [21 ] on PNSs to consider differential rotation and uni¬ 
versality. One could also study whether the approximate 
universality found in the new four-hair relations could 
possibly be improved through a different choice of nor¬ 
malization, which was recently found for the three-hair 
relations in the absence of differential rotation [36]. 
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Appendix A: Stellar Shape Deformation for 
Differentially Rotating Stars 

We here extend the analysis of Sec. 2.4.4 in [75] to 
derive the stellar shape of differentially rotating stars to 
leading (Newtonian) order in a weak-field expansion. We 












11 


work in a perturbative differential rotation and a slow- 
rotation approximation, and thus, we expand the shape 
about a non-rotating background. The background sat¬ 
isfies a mass conservation law and a hydrostatic equilib¬ 
rium equation: 


m(r) = f p(r)r 2 dr , ^ = - p(r ) m ^ , (Al) 

Jo dr r 

where R is the stellar radius in the non-rotating config¬ 
uration. 

Let us next look at the perturbed equations. Focusing 
on a generic, non-axisymmetric perturbation, the per¬ 
turbed mass conservation law and hydrostatic equilib¬ 
rium equation are given by 

6 p = —£ j djp f ^djp- -djSp + dj(5U+ V) = 0, (A2) 

(r p 

where Sp, 5p and SU are the perturbations to the pres¬ 
sure, the density and the gravitational potential respec¬ 
tively. We have also introduced as the displacement 
of a fluid element in the perturbed configuration relative 
to the background and V as the potential that drives the 
perturbation (namely, the centrifugal potential). 

Let us now decompose the perturbed quantities in 
spherical harmonics: 


€ = Y, rf em (r)Y im (0, 4 >), (A3) 

tm 

Sp = '^2p(m(r)Y em (d, </>), (A4) 

tm 

&U = Y d Utm(r)Y trn {p,<l>), (A5) 

tm 

V = Y,Ve m (r)Y ern (6,<f>). (A6) 

tm 


Substituting the above decomposition into the perturbed 
mass conservation law in Eq. (A2) and using the EoS 
and the background hydrostatic equilibrium equation in 
Eq. (Al), one finds 


pm 

Ptm = 

r 


ftm • 


(A7) 


and thus, we need to derive a relation between f £m and 
Vi m . From the perturbed Poisson equation in the exterior 
region, one finds 


Ue m (r) 


47r M £m 

2 £ + 1 r*+! ’ 


(All) 


where M £m correspond to the stellar mass multipole mo¬ 
ments. For an axisymmetric configuration, M £ q agrees 
with M £ in Eq. (1) modulo a constant factor. Equa¬ 
tion (A9) at r = R then becomes 


R fem{R) 


47r M lm 
21 + 1 RC+ 1 


+ Yim{.R) i 


(A12) 


and a similar equation holds for the first derivative of 
Eq. (1). Eliminating M £m from these two equations, one 
finds 

hm{R) = yyy y [(* + WM + RVLm , 

( Ai3) 

where we introduced the surficial Love number h £ , de¬ 
fined by 


21 + 1 

f + vt{R) 


(A14) 


which depends on the EoS, and the Radau function 


m(r) 


r fL 

i im 


(A15) 


Let us now focus on an axisymmetric rotational pertur¬ 
bation and consider the m = 0 mode only. Expanding in 
X = Si/Mq <C 1 and in |-y| 1, the centrifugal potential 

is given by 


V(r,0) 


^fi(r, 0) 2 r 2 sin 2 9 + 0 (\ 4 ) , 

\ny [i +pm 


8 n 2 y 
15 7 ^T 


10 ^ 

l-yW + ^W 


+ 0 ( X \j 2 ), 


(A16) 


Similarly, from the angular component of Eq. (A2), one 
finds 


where recall that p, = cos 9 and we used Eq. (3) for 
f l(r,9). Comparing Eqs. (A6) and (A16), one can read 
off Wo for the rotational perturbation: 


Ptm — p(Utm H - Vtm) • 

(A8) 


Combining Eqs. (A7) and (A8), one finds 


Wo — 

7TI 

ftm = tm “1“ Vtm • 

r 

(A9) 

Wo = 

The goal of this appendix is to derive the stellar shape 
deformation, 

Wo = 

5R = Y J Rhm{R)Y £m {0A), 

tm 

(A10) 



2 ^/tt 
3 

2 V / 5n : 

15 

16 s/tt f l c 2 r 4 
105 7 ai 2 


SiA 2 ( 1 -§7p)+0(X,7 2 ). (A17) 

SicV (i-y7 


A ,~,2\ 


(A18) 

(A19) 
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From these expressions, one can obtain the perturbations 
and from that the stellar shape deformation us¬ 
ing Eqs. (A10) and (A13). Re-expressing R and fl 2 in 
terms of ai and e 2 , one finds 


SR = or <{ 1 — - [1 + 2 P 2 (/i)] e 2 


16 

H- 7 

15 ' 


p 2(v) - jrrr p ^) 

14 h 2 


2 }+0(e 4 ,' 


(A20) 


The uniform rotation part in the above equation corre¬ 
sponds to the spheroidal shape deformation, while the 
0 ( 7 ) part is a correction due to differential rotation. 

Once the stellar shape deformation has been found, one 
can read off the Legendre coefficients a 2 j- Comparing 
Eq. (A20) to Eq. (9), one finds 


Since he enters only at 0 ( 7 ), we need to evaluate it, and 
thus the Radau function, for a uniformly rotating star. 
The Radau function satisfies the Radau equation, which 
for uniform rotation is 


3 

rq'e + r]e(rie - 1) + 87 t— (rj e + 1) - 7(7 + 1) = 0, (A22) 
m 

with the boundary condition 7 ^( 0 ) =1—2. For an n = 
0 polytropic EoS, and taking into account the density 
discontinuity at the surface, one finds 


,(n= 0 ) 27+1 

h > = W^T)' 


(A23) 


Substituting this into Eq. (A21), one obtains ae for an 
n = 0 polytrope: 


1 fi 

a 0 = 0 (e 4 ), a 2 = —e 2 + 0 (e 4 ), 

15 

94 

a 4 = - — e 2 + 0(e 4 ). (A24) 


a 0 = O (e 4 ) , a 2 = ^e 2 + O (e 4 ) 


CH 4 = 


8 /14 

35 ho 


e 2 + 0 (e 4 ) . 


(A21) 


Appendix B: Detailed Expressions for the Multipole 
Moments for Differentially Rotating Stars 


We here present detailed and generic expressions for the mass and mass-current multipole moments with differential 
rotation within the perturbative differential rotation and weak-field approximations. Normalizing Eqs. (16) and (17), 
we find 


M 2 e + 2 = 


(-i ) t+1 3 2 *+ 2 £f +4 IoT \W{Hl)\ n+1 Rn,U + 2 

4 ^+! f> 2 *+ 2 af f2 6 I ?+ 2 


+ 4 (e + i) SIl ’ 5 ~ SI] _ 

3 + j 5/1,3 i\Rn, 


T - 1 

3 1 2i+2,3 

1,5 R-n, 4 


77; 


2^+2 


1 + 7| 3(27 + 1) ^ - 10(7 + 1) ^ + (27 + 5 ) Iu+2,2 


*0,3 


5 / 1,3 


o(Y 


/2t+2,3 

(Bl) 


S 2 t +1 — 


(-l)'3 2 ' +1 2f+Hf I$‘ 3 |tf'(|,)| 2, R„, 2 , +2 
4< Si r it+3 a? SIS' J/ 2 7 +u 


1+7 


67 


ki^ 5{2i+1) s Ai + {2i+5) s J^±hi 

m,3 o7l,3 0121+1,3 


2(21 + 1) 


5/l,5 — 5 /i, 5 77„,4 7 + 1 5/2^-(-i,5 — 572^+1,5 77 n ,2/+4 


357 13 Ci77„,2 47 + 3 blii+i,?, C 2 77„,2^+2 / _ 


0 ( 7 2 


(B2) 


These multipole moments are the basis for the universal relations that lead to the three- and four-hair relations in 

the presence of differential rotation. The first and second universal relations of Eqs. (21) and (22) are corrected due 

(1 2 ) 

to differential rotation through the coefficients oc n \ , which are explicitly 


a (1) - 

a nl 


! 

2 Vl^' 


557^2 + 157+f + (-l) f+1 (27 + 3) (27 + 5) 


ht+2,2 27+ 1 572^-1-1,2 

e 2e + 2 47 + 3 e 2f 


( 2 ) 

al 0 = 


4 (_l)f-K(£ + i)(27 + 1) (2e 2 + 7) - 

(27 + 5)(47 + 3)e 2 n,t ’ 

.7 + 1 


2 Vl zr e 2 


-15^’ 


n, 0 


5ii,2 — 3 / 0,2 “I - ( — + 3)(2/ + 5) 


l2t+2 ,2 7 + 1 27 + 1 5 / 2 ^ 71,2 


0 2£+2 


7 47 + 3 e : 


2£ 


15 7 


+ 1 4(—1) £+1 (7 + 1)(27 + 1) (2e 2 + 7) 7 + 1 - 

C n ,Q H-„ .—_ w . „ . n-7—bn 


(27 + 5) (47 + 3)e 2 


(B3) 


(B4) 
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The three- and four-hair relations of Eqs. (24) and (27) in the presence of differential rotation are corrected by the 
coefficients and a^’ M,S \ which are explicitly 


a 


(3 ,M) 


2vT^ 


3£/o,2 + M>{£ + 1) —i—b ( 

e 2 


-lY +1 {2£ + 3){2£ + 5) ! ^g^ 


(B5) 


x n,tL 


2 v / )T^ 


3£/ 0j2 + 15*% - 561 li2 + (-1)*(2* + 3)(2* + 5)§^^tM 
e- 4e + 3 e 2t 
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